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1. Introduction 

At low temperatures and densities the quarks and the gluons are "confined" in colorless states 
such as the nucleons and the mesons. In other words, there are no free quarks or gluons detected 
experimentally. The question of what happens to ordinary hadronic matter at extreme conditions 
(very high temperatures or densities) has been recently answered by the heavy-ion collision (HIC) 
experiments at RHIC and LHC. At these experiments, heavy nuclei accelerated close to the speed 
of light collide to form a new state of matter called the quark-gluon plasma (QGP). Prior to its 
discovery, the existence of the QGP had been speculated about for more than 30 years [|] on the 
basis of the expectation of asymptotic freedom of QCD [Q] at high energies. The properties of the 
QGP are important for the study of a variety of physical phenomena. For example, up to several 
microseconds after the big bang all the hadronic matter in the Universe was in a form of QGP. 
Hence, the equation of state (EOS) of QGP plays a role in the expansion of the early Universe 
and may influence the subsequent structure formation. The QGP also may exist in the cores of 
neutron stars (at very high baryon density) and may even form quark stars. And lastly, for HIC 
experiments the EOS of QGP is an essential input to the hydrodynamic equations used to simulate 
the expansion of the created plasma. The initial temperatures occurring at the HIC are not more 
than five times higher than the transition temperature T c between ordinary matter and quark-gluon 
plasma. In this temperature range the quark-gluon plasma is in a nonperturbative regime and only 
a nonperturbative tool, such as lattice QCD, seems to provide an adequate description. 

Lattice QCD is simply the discretized version of QCD on a space-time lattice of dimensions 
Nf x N t . The temperature on the lattice is defined as T = 1/ (aN t ), where a is the lattice spacing. To 
change the temperature on the lattice one can fix N t and vary a. This type of temperature variation 
preferably should occur on lines of constant physics (LCP). A different strategy would be to fix 
the lattice spacing and change T by varying N t . This way staying on an LCP is automatically 
guaranteed. The continuum limit at a particular T is taken as N t — > °° (i.e., a — > 0). Generally, 
to insure that a lattice calculation will not suffer from large finite volume effects, the following 
inequality should be fulfilled: 

a<M / ; 1 <aN s , (1.1) 

where M/, is the mass of some characteristic light hadron. The temperature on the lattice can be 
considered high when TM^ 1 ~ 0(1), which means that the above inequality in this case is modified 
to: 

l/N t 4Zl<.N f /N t . (1.2) 

Current simulations are done with N t = 6 — 16 and aspect ratios N s /N t = 3 — 4. The discretization 
effects at finite a will strongly depend on the choice of lattice action. There are two main types of 
fermion lattice actions in use: staggered types [Q] and Wilson types [J3j] (including chiral formu- 
lations [Q]), with the former being the choice with the largest investment of computing resources 
thus far at high T. 

2. Staggered actions and discretization effects 

The most advanced nonzero temperature calculations so far have been done with staggered 
fermions. This type of fermion discretization is one of the most computationally affordable and 
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Figure 1: (Left) The free energy normalized to the Stefan-Boltzmann value vs. the square of the inverse of 
temporal extent in the free theory case. (Right) The square of the RMS pion mass vs. the lattice spacing 
Jzo|], It is assumed that the Goldstone pion mass is tuned to its physical value. The vertical lines denote the 
approximate lattice spacings at which the transition happens for different A/f's. 



thus it still remains as the most popular choice for the expensive high-temperature calculations. 
Another advantage of the staggered fermions is that they retain af/(l)xf/(l) symmetry (for one 
staggered flavor), which is a remnant of the full chiral symmetry group. This remaining symmetry 
protects the chiral limit from an additive renormalization in the quark mass. The main drawback 
of the staggered type of discretization is that for a single staggered flavor there are 4 "tastes" of 
Dirac fermions appearing in the continuum limit. To deal with this taste proliferation, the usual 
strategy is to work with the fourth root of the fermion determinant instead of the full one in the 
Monte Carlo generation of gauge ensembles. The hope is that this procedure effectively reduces 
the number of tastes to one per flavor, although this has not been proven rigorously. Still, there 
is a growing amount of numerical ^ and theoretical [f7|] evidence that the fourth root procedure 
gives the correct number of flavors in the continuum limit, provided it is taken before the chiral 
limit. At finite lattice spacing the taste symmetry is broken, which means that the staggered hadron 
multiplets are nondegenerate. For example, there are 16 staggered pions and the fact that they have 
nondegenerate masses can distort the physics at not-too-high temperatures below the crossover, 
where they are expected to dominate it. 



There are a number of improved staggered actions in use: the p4 \f$ft, asqtad [Bp, stout QTQp and 
the HISQ [11] action, which deal with the discretization effects in somewhat different ways. The 
pA action has an improved quark dispersion relation to 0(p 4 ), where p is the quark momentum. 
The asqtad action also has this property but furthermore, it has all tree-level lattice artifacts to 
0(a 2 ) removed, and its coefficients are tadpole-improved. The stout action has the taste symmetry 
violation reduced through unitary link smearing. And finally, the HISQ action combines all of the 
above improvements. Comparisons of the sizes of the lattice artifacts among the above actions are 
shown in Fig. [T[ In the left panel the approach to the continuum limit of the free energy in the 
free fermion case is presented. The stout action has a slower approach to the Stefan-Boltzmann 
continuum value than the rest. The closer the value of the free energy from a particular action is to 
the continuum the smaller are the UV cutoff effects for this action in the high-temperature regime 
above the phase transition. However, in the regime of low temperature or around the transition, the 
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Mf = 2 Pure 




Figure 2: (Left) The QCD phase diagram at zero chemical potential. (Right) Fits with 0(N) scaling mag- 
netic equations to the chiral order parameter p0| ] . 



taste symmetry violation would seem to play a more important role in determining the discretization 
effects. If we choose the root-mean-square (RMS) pion mass as a measure of the taste symmetry 
violation defined as 

m* MS = J -j^r (m\ + m\ yi + 3m| ^ + 2>m\ y . + 3m| % + 3m| + m\ + mf) , (2.1) 

where the sum is over the squares of the masses of the members of the staggered pion multiplet, 
then the HISQ action appears to have the smallest taste symmetry violations of all at any a's. On 
the other hand, the stout action outperforms the asqtad action for coarser lattices. 



3. QCD at zero density 

The widely accepted conjecture for the QCD phase diagram in the plane of the strange quark 
mass m s vs. the light quark masses m u ^ is shown in Fig || (left panel). At a range of low or high 
quark masses (shown as green areas in the lower-left and upper-right corners), the expectation is 
that the phase transition at high T is of 1st order. In the first case it should have the character of 
chiral symmetry restoration and in the second — of a deconfining type. Both of these areas are 
surrounded by critical lines of 2nd order transition. At zero m u ^ and sufficiently large m s , the 2nd 
order critical line changes from Z(2)-type to 0(4) -type at the tricritical point mf c . The current 



consensus is that the physical point is somewhere in the vast area of the crossover region []12|]. It is 
important to study the phase diagram quantitatively, especially when locating the physical point in 
relation to the low-mass region of 1st order transition. If the physical point lies low, when taking the 
chiral limit one encounters a Z(2) 2nd order phase transition at nonzero mass. On the other hand 
when it lies higher, the 2nd order transition is expected to be of 0(4)-type. Establishing which 
particular behavior is true determines the type of scaling functions one should use to describe data 
at low quark masses. In a recent study [ p3| ] with 2+1 flavors of stout fermions, it was found that 
the 2nd order critical behavior starts at quark masses <12% of the physical quark masses. This 
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Figure 3: (Left) Comparison of asqtad data for X m ,l with the differentiated scaling O(N) curves ||2C|]. (Right) 
Comparison between the HotQCD and WB data for the subtracted chiral condensate [Effl], The continuum 
and physical quark limit for the WB data is denoted by pink triangles. The same for the HotQCD data is 
denoted with black diamonds. 



suggests that the physical point is relatively far from the region of the 1st order transition and the 



chiral phase transition is probably of 0(4) -type. This year, a new study [14] with 3 light degenerate 



flavors of HISQ fermions reached similar conclusions — the critical quark masses are found to be 



< 10% of the physical light quark mass. Another new study Q15| ] with 3 degenerate flavors confirms 
the Z(2)-type of scaling behavior as expected from the qualitative picture of the phase diagram. 
The upper-right corner of the phase diagram has also been explored numerically. The WHOT- 



QCD collaboration has done a recent study Qlq ] of the heavy-quark region starting from pure SU (3) 
configurations (i.e., infinite quark masses) and using reweighting combined with hopping parameter 
expansion to reduce the quark masses and map the critical line. This year they extended this work 



by adding nonzero chemical potential [17]. 

3.1 QCD and universal scaling 

As described in the previous subsection, there is an increasing amount of evidence that in the 
chiral limit of 2+1 flavor QCD the phase transition is of 2nd order in the 0(4) universality class. 
Close to the chiral transition (in temperature and quark masses) the free energy can be written as a 
sum of a singular part f s and a regular part f r : 

L(r,m u m s ) = h w ^f s (z)+f r (T,H,m s ), (3.1) 

where z = t/h 1 ^^ 8 ^ is a scaling variable and /3 and 8 are critical exponents of the 0(N) universality 
class. The reduced temperature t and the symmetry breaking field h are defined below: 



1 t 



l(l^ +LKq (tk)\ Kls ^) , h = L H , H = ^ ( 3.2) 
'o V r ° q =i, s ) ls T T ) h m ; 



where in the general case t depends on the light and heavy quark chemical potentials \i[ and \L S (in 
the 2+ 1 flavor case) and T c is the true chiral critical temperature. The dynamics around the phase 
transition is governed by the scale invariant singular part f s , which means that to extract the critical 
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temperature it is best to choose observables sensitive to it (the regular part f r gives rise to scaling 
violations in these observables). Such an observable is of course the chiral order parameter, which 
has the critical behavior (i/n//)/ ~ m 1 / 5 at fixed z. Another observable sensitive to f s is the chiral 
susceptibility %m,l = d{\jfyf)i/dmi with the critical behavior Xm,l ~ m i ■ The % m j is divergent 
in the chiral limit (8 > 1 for 0(N) scaling) and the position of its peak is often used to define 
the critical (or crossover) temperature. An observable also sensitive to f s , but not widely used, 

is the mixed susceptibility %tj = —T d 2 \nZ / {V dmidt), which has a divergent critical behavior: 

(B-l)/(BS) 

To define the critical temperature, sometimes so called "deconfinement type" observables are 
used, although their relation to the 0(N) critical behavior is less clear or absent. One such ob- 
servable is the quark number susceptibility Xq=l,s/T = d 2 lnZ/(VT 3 d(iu. q /T) ) which quantifies 
the light or strange quark number fluctuations. The Xq=ls are expected to be zero at T = 0, where 
the quark flavors are confined in colorless states, and to increase as the temperature grows and the 
quark flavor is "freed". The quark number susceptibilities are also important for the study of HIC 
since they are connected to event-by-event fluctuations of various quantum numbers such as charge, 
strangeness and baryon number. For more details on how Xq=l,s can be connected with various ex- 
perimental observables see Ref. [18]. The quark number susceptibilities are somewhat sensitive to 
f s but they do not become divergent in the chiral limit. For example the critical behavior of the 
temperature derivative of Xq i s: 



d Xq ^ A 



T — Tx 



(3.3) 



where a < for 0(N); hence, it is not divergent, and in fact, it is dominated by the regular part c r . 
Thus 0(N) scaling fits to the quark number susceptibility are not very reliable for extracting the 
critical temperature. 

Another often used deconfinement-type observable is the Polyakov loop which is related to 
the free energy of a static quark: (L) ~ exp(— Fq/T). In the limit of infinite quark masses it is 
governed by the singular part of the partition function of the pure gauge theory, and it is an order 
parameter for deconfinement. At finite quark masses it is not, and its use as a mean to define and 
extract T c has a less firm theoretical justification. Hence in this review I will concentrate on results 
for the critical temperature obtained through only the chiral-type observables, which are expected 
to be sensitive to the universal 0(N) scaling. 

The extent to which universal scaling is applicable has been studied with both staggered 
[ ]19i 20] and Wilson [|T]] fermions. Generally it appears that the 0(N) scaling gives an adequate 
description of the data. At finite lattice spacing the scaling for rooted staggered fermions should 
be 0(2) -type rather than 0(4), but in practice the critical exponents for both types are similar 



enough that numerically they cannot be well distinguished gOJ . To illustrate this statement, in 
Fig. |2| (right) scaling fits from Ref. J20| ] to the renormalized order parameter Mb = m s (\fir\ff)i/T 4 = 
h 1 / 8 ' fc{z) are shown. Both 0(2) and 0(4) magnetic equation of state fits appear to perform simi- 
larly well. The resulting extrapolation to physical quark masses is also shown. If scaling is applica- 
ble, then by differentiating these scaling fit curves with respect to the light quark mass, one should 
be able to predict the position of the peak of Xmj- In Fig. |3] (left) a comparison between the Xm,i 
data and the differentiated scaling curves is shown. Indeed, the predicted and measured positions 
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of the peaks line-up satisfactorily. The scaling curves do not include regular parts and this is the 
most likely reason for the discrepancy between the %m,l and the curves at lower T. 

3.2 Critical temperature 

Since at physical quark masses there is no true phase transition to QGP at high T, the definition 
of a critical temperature T c is to some extent a matter of choice. However, a determination of T c 
based on observables sensitive to the singular part of the free energy f s is on a firmer theoretical 
footing and has a clearer interpretation. For this reason all of the results for T c quoted below are 
extracted from the peak in Xm,l- 

Historically there appeared to be some discrepancies between T c values calculated for 2+1 
flavors with different types of improved staggered fermions. In 2005 the MILC collaboration re- 



ported [g2J T c = 169(12) (4) MeV from a calculation with asqtad fermions. In 2006 the BNL- 



RBC-Bielefeld collaboration published [gj the value 192(7) (4) MeV, using the pA action. The 



Wuppertal-Budapest (WB) group in the period 2006-2010 consistently obtained [ |24j ] lower values 
with stout fermions, the latest of which is 147 (2) (3) MeV. All of the above numbers are extrap- 
olated to the continuum limit. The MILC result is also extrapolated to the chiral limit, while the 
other two values are at the physical quark mass limit. The latest number from the HotQCD collab- 
oration, obtained from calculations J2CJ| ] with both the asqtad and the HISQ action is 154(9) MeV 
(extrapolated to the continuum and physical quark mass limit). This result is compatible with the 
WB value. To illustrate this statement in Fig. ^] (right) the continuum and physical quark mass 
extrapolations of the WB and HotQCD data for the subtracted chiral condensate 

A/, s = ((yy)i.t - (yy)s,t—)/((VY)i.o ~ (YY)s,o—) (3.4) 

are shown [p0|]. The stout and HISQ data have practically identical extrapolations, suggesting the 
same T c . As for the higher value of the BNL-RBC-Bielefeld result, it is probably because the 
calculation was done at TV, = 4 and 6 only, which means relatively coarse lattices were used and the 
discretization effects proved to be significant. 

All of the above results for T c have been obtained in the limit of "infinite" volume approximated 
through enforcing (anti)periodic boundary conditions (PBC) on the lattice. At HIC experiments 
the volume of the plasma created is only 5-10 fm 3 [25]. The effect of such a small volume on the 



phase transition may be significant. In 2007 there was an attempt [26] to estimate it on the lattice 
by calculating in pure SU (3) the effects of cold boundary conditions (CBC) on T c (meaning that the 
ends of the lattice were left "free"). This setup is supposed to imitate a hot plasma immersed in a 
cold (zero temperature) exterior. The results are shown in Fig. || (left) where the CBC are compared 
to the periodic ones at different spatial lattice sizes L s . It seems that the corrections to T c due to the 
finite volume effects can reach up to 30 MeV. However, the CBC are not too easy to interpret, since 
the exterior temperature cannot be determined straightforwardly. For SU (3) in the scaling regime, 
the temperature can be set through the RG equation, but the cold exterior has essentially zero 
coupling which is out of the range of its applicability. To address this difficulty this year an extended 



study [27] was presented with a setup where both the interior and the exterior are kept in the 
scaling regime of SU (3), through the so called double-layered-torus (DLT) boundary conditions. 
The exterior has temperature of around 158 MeV, while the interior is at T c (as determined from 
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Fits of ff translated into T°(MeV) (T OJ =158.15MeV for DLT) 




Figure 4: (Left) The effects on T c from different lattice boundary conditions p7|]. (Right) The dependence 
of T c on the strength of the magnetic field eB obtained from the peak of the average of the renormalized 
chiral condensates %' ud [p4|]. 



the Polyakov loops). Even with such a hot exterior by comparison with the HIC conditions, the 
T c with the DLT is still 6-12 MeV higher than in the infinite volume case. It will be interesting to 
know what the finite-volume correction would be with dynamical fermions added. 

Another factor that may influence the value of T c is the presence of strong external magnetic 
field. At HIC the magnetic field created in noncentral collisions is of O(10 14 ~ 15 T) [fZ8|], a value 
much larger than a field of a typical magnetar O(10 10-11 T) [ j29| ] , and not too far from the fields 
created at the beginning of the Universe — O(10 16 T) pOfl. Strong external magnetic field may 
have a pronounced effect on the QCD phase transition, not only by modifying T c but also by 
changing its order, "separating" the chiral restoration from the deconfinement and inducing CP- 



violation through the Chiral Magnetic Effect (CME) [|31[]. Some effective models [ ]32[ ] predict that 
T c will increase with stronger external magnetic fields. Others find the opposite effect [33]. It is 
important to settle these discrepancies with a first principle lattice calculation. This year new results 
[34] for the magnetic field effect on T c with stout fermions with physical masses 1 became available 



and are shown in Fig. ^ (right). The calculation was done at N t = 6, 8 and 10 and the continuum 
limit was taken. The results show a decrease of T c when the magnetic field eB increases, which 
is the opposite of what was found in Ref. [|35|]. However, the latter calculation was done at heavy 



quark masses (which the authors of Ref. [ ]34[ ] claim as the most important factor for the difference), 
coarse lattices and with standard staggered action, all of which may contribute to the discrepancy 
with this year's more advanced result. In both of these lattice calculations the setup is such that the 
fermions interact with the external magnetic field but do not interact electromagnetically among 
themselves. This "quenched EM" approximation introduces an unknown error in the final results. 
In the new study [|34j] also no evidence was found for the transition becoming a 1st order in the 
range of the magnetic fields available (up to around 1 GeV 2 ). 



'When it comes to staggered fermions, "physical" quark masses are defined as the ones which give the lightest of the 
staggered pions and kaons their physical masses. It can be argued that since the rest of the staggered multiplet members 
are heavier, the effective quark mass is actually larger than the physical even in this case. 
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Continuum extrapolation of s/T 
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Figure 5: (Left) The trace anomaly calculated with different staggered action and compared to the HRG 
model at lower T [fto|], (Right) The entropy density calculated at three values of T with the new method 
from Ref. 



3.3 The equation of state 

The QCD equation of state (EOS) gives the temperature dependence of energy density, pres- 
sure and other bulk thermodynamic quantities of a quark-gluon system. It is an essential input to the 
hydrodynamic equations used to model the expansion of the plasma created at HIC. In recent years 
the demand for the EOS calculated from first principles has grown [|^] as the HIC experiments 
accumulate high statistics on spectrum data. Lattice determinations of the EOS are computation- 
ally expensive and historically the most resources have been invested in the integral method with 
variable scales. The integral method is based on the thermodynamic identities for energy density 
£, pressure p and the trace anomaly 0^: 



eV 



dlnZ 



d(l/T) 



P_ 

T 



dlnZ 



dV 



InZ 



3p 



Td\nZ 
V d\na 



(3.5) 



From the above it follows that p and e can be obtained by integrating 0^ over the lattice spacing 



(hence the method's name). In the most widely used version of this method [37], the lattice spacing 
is varied along an LCP in order to change T while N t is kept fixed. This means that at lower 
temperatures the lattices are coarser and the discretization effects larger. In the fixed-scale version 



of the integral method ||38|], where a is kept constant and N t varies (and the integration is over N t ), 
the situation is reversed — the larger discretization effects will be at high temperatures. There is 
another fixed-scale method which does not require integration, called "the operator/differentiation 
method" [p9|]. The pressure and energy density in this case are calculated directly from the first 



and second expressions in Eq. ( p75| ) (i.e., not through the determination of W first). However, 
the method demands the determination of anisotropy coefficients in addition to the beta- and mass- 
renormalization functions required for the integral methods. All of the methods described require 
a zero temperature subtraction to eliminate the UV divergences in the EOS. 

In Fig. |5] (left) the values for the trace anomaly are compared for several staggered actions and 
at different A^'s. The method of calculations is the variable scale one. The data for the asqtad, p4 



and the HISQ action is from the HotQCD collaboration [4C]. It is calculated on an LCP, where 



m s is physical and mi/m s = 0.05. The stout data is from the WB group pi] ] and is obtained with 
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Figure 6: (Left) The pressure obtained with the heavy-quark quenched approximation in the 2+1+1 flavor 
case and compared with the 2+1 (dynamical) flavor result [55]]. (Right) The pressure calculated with 2+1+1 
dynamical flavors compared with the 2+1 flavors [E6|] . 



all quark masses set to the physical ones and at N t = 8 (they also have data for N t = 10 and 12 
|gl|, which agree with the N t = 8 case). At lower temperatures all the results are compared 
with the hadron resonance gas (HRG) values, which is expected to be a reasonable description of 
the system in this regime [42]. The data at larger N t (finer lattices) start to approach the HRG 
result at low T and the small difference which remains between the WB and HotQCD data there 
may be attributed to the difference in the quark masses. The more striking difference in the region 
above the crossover is still without explanation. However it is probably not a result of quark-mass 
differences (or deviations from the LCP), since these should play a minor role at high T. It is 
possible that additional HISQ runs at N t = 12 may help to explain this discrepancy. 

For HIC at RHIC the EOS at 2+1 flavors is considered most relevant, since the time scales are 
probably too short for the thermalization of the charm quark [43]. At LHC the quark-gluon plasma 
is hotter and lasts longer and thus it is possible that the charm quark thermalizes in this case. This 
implies that its contribution to the EOS of the plasma created at LHC may not be negligible. Pre- 
vious heavy-quark quenched calculations [|4|, ^3]] (i.e., where the sea charm loops are neglected) 
attempted to calculate the effects of the charm quark (see Fig. || (left) for pressure results). But 
the question remained whether the quenched approximation used for the charm quark introduces a 
substantial systematic error. A new study J+6| ] with stout fermions at N t = 6,8 and 10, where all 
the u,d,s and c flavors are dynamical, points toward an answer to this problem. In Fig. |6] (right) the 
results for the pressure for 2+1+1 dynamical flavors is presented. By comparison with the heavy- 
quark quenched result from Ref. [45] in the left panel of the same figure, the contribution of the 
charm quark is smaller throughout and becomes significant only at relatively larger temperatures 
(T > 350 MeV). The authors also show that dynamical c-quark effect in the pressure is very close 
to its perturbative estimate [47]. They attribute the difference between the quenched and dynam- 
ical cases to a shift in the LCP which results from the c-quark being part of the sea. There still 
remains the question of the heavy-quark discretization effects and how well the charm quarks are 
represented on the lattices used in Ref. [46] with stout fermions. 

The EOS is also studied with Wilson-type fermions; Fig. [7| shows two of the most recent 
results. The WHOT-QCD collaboration used the integral method with fixed scales and the data 
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Figure 7: (Left)WHOT-QCD EOS result iQ for 2+1 flavors of nonperturbatively improved Wilson 
fermions. The lattice scale is fixed to a « 0.07 fm and N t is varied between 4 and 16 and m K /m p « 0.63. 
(Right) New result J48| ] for the trace anomaly calculated with 2 flavors of mtm Wilson fermions at N t = 12. 
The lattice scale is varied 0.04 — 0.1 fm. The outliers are probably due to finite volume effects and zero- 
temperature interpolations. 



from Ref. [38] is shown in the left panel of this figure. In the right panel of Fig. [7] a new result [48] 
with mtm Wilson fermions is presented. In this case the method with variable scales is employed. 
More details of these two calculations are given in the caption of Fig. [7[ Both of these studies are 
at heavy quark masses and have large statistical errors, and for the time-being, they are not yet 
as sophisticated as the staggered calculations. However, Wilson studies are essential for verifying 
the current staggered results with an alternative action which does not suffer from taste symmetry 
violations. This year such a staggered- Wilson comparison was reported for some thermodynamic 



quantities [49] (namely: (yy), % s and the Polyakov loop), which show a very good consistency 
between the two different types of fermion action. 

A qualitatively new method for calculating bulk thermodynamic quantities was recently intro- 
duced in Ref. [ p0| ] . It is based on a lattice determination of the generation function of the cumulants 
of the momentum distribution K(/5,lJ,,z)- The cumulants of K(p,n,z), where z is a 3 component 
vector, defined as 

d 2, n d 2n 2 d 2n 3 K(j3,H,z) 



^2ni,2n 2 ,2n 3 — ( — 1) 



dz?" 1 dzl n2 dz]" 3 



\z=0 



are related to the entropy density s and the heat capacity c v in the following manner: 

^2,2,0 3^2,0,0 



K 



2,0,0 



(3.6) 



(3.7) 



J A J2 v ' 

The above means that the knowledge of K(p,/J.,z) would give access to the determination of the 
EOS, since by integrating s over the temperature one can obtain the pressure using the identity 
s = Tdp/dT. To calculate K(p,pi,z) the following result from thermal quantum field theory is 
employed: 



where Z(j8,/i,z) = Tr{e~P( H ~ tiN ^~' pz } is a partition function where all states with momentum p 
acquire a phase e' p z . On the lattice this partition function can be represented as a Euclidean path 
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integral with shifted b.c. for the fields in the temporal direction: 0(/,Jt) = ±0(O,x + z). According 
to the method in Ref. [50] applied at pL = 0, to evaluate the partition function ratio in Eq. (3.8) a 
system of N interpolating ensembles with Z(j8,r,) is created, each with the following action: 

s(u,n) = n s(u) + (\- ri )s(u z ), (3.9) 

where S(U ) is an action with PBC, S(U Z ) is one with shifted b.c. and r,- = i/N, i = 0, . . . , N. To have 
the discretization effects under control the number ,/V should be 0(N^). Then K is approximated 
as: 

Wz,a) = -m^f « -£h* = - Iln(e^-)-^)), +1 . (3.10) 
The derivative of K(p,/j.,z) to determine the entropy is taken as the following limit: 

jW s M =h y, (3 . U) 

T 3 T 5 a^O \z\ T 5 N^ 

The deviation from the continuum at finite lattice spacing with the above procedure is 0(a 2 ). The 
authors of Ref. [50] have applied the method to the pure St/ (3) case using a plaquette action. 
In Fig. H (right) the results for s on the lattice spacing are shown for three temperatures. The 
approach to the continuum at the lower two T's seems rather mild. At the highest available T the 
discretization effects seem more pronounced, probably due to the lattice spacings being somewhat 
coarse for T = 9.2T C and the UV cutoff being more pronounced in this case. The computational 
cost for this new method could be as high as 0(a~ n ) for the pure gauge theory depending on 
the algorithm for estimating the ratio in Eq. (|3lj|). The method does not require zero temperature 
subtractions to eliminate the UV divergences in thermodynamic quantities, but the overall cost is 
likely not going to be lower than the one for the more traditional methods. 



4. QCD at nonzero baryon density 

Studying QCD at nonzero baryon density (or nonzero chemical potential /i) is important for a 
variety of physical phenomena such as supernova explosions, neutron star formation, the existence 
of quark stars and degenerate matter. The properties of hot hadronic matter at nonzero baryon 
density are also explored experimentally at lower-energy HIC. Thus a quantitative determination of 
the QCD phase diagram in the T— baryon density plane is of both great theoretical and experimental 
interest. The current conjecture for the diagram is shown in Fig. || (left). At high T and low density 
it is accepted that the transition is a crossover. It is expected that as the baryon density increases the 
crossover region moves to lower T until a 2nd order critical end point (CEP) is reached. At even 
higher densities the transition is supposed to become of 1st order. On the diagram the first order 
region is denoted by a band since for a 1st order transition the two phases coexist at a (narrow) 
range of baryon densities at the critical temperature. This first order band shrinks to a line when the 
diagram is presented in the T — fi plane instead (i.e., the diagram in Fig. |8| (left) is from a canonical 
ensemble point of view, instead of the grand canonical one in the T — pL plane). It is also important 
to establish quantitatively the relative position of the experimental freezeout curve [ |5l| ] with respect 
to the critical one, since if the two are far apart the measured signals at HIC experiments for the 1st 
order phase transition may become "washed out" and more difficult to interpret. 
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Net Baryon Density 



Figure 8: (Left) The QCD phase diagram at nonzero baryon density. (Right) Estimate of the phase boundary 
at T = 190 MeV (which is slightly higher than Te) as a function of the order of the expansion from r^, r% 
and Pade approximations of order p/q J68[|. The dotted line is the full model solution. 



4.1 Nonzero chemical potential on the lattice 

Introducing nonzero chemical potential on the lattice presents a significant challenge — the 
straightforward Monte Carlo (MC) simulations, where the determinant of the fermion matrix detM 
is interpreted as a probability density, become unfeasible. The reason for this is that detM(/x ^ 0) 
becomes complex: 

(detM(ju))* = detM(-/i) ^detM(/i). (4.1) 

This constrains the MC simulations to pL being zero, or pure imaginary, or to the case of pair(s) of 
degenerate flavors with opposite real pL (nonzero isospin chemical potential) since in those cases 
detM remains real. There are a host of methods which attempt to avoid the problem of the complex 
determinant (also referred as "the sign problem") at jj. / 0. Examples are: the Taylor expansion 



method [|52|], reweighting techniques [p3Q, analytic continuation from imaginary \x [[541], density of 
states methods [ |55| ] and canonical ensemble approaches [ |56| ] . All of these methods find a way to 
circumvent directly simulating at real nonzero /x. In the next subsection I discuss in more detail 
the Taylor expansion method only. For recent reviews of the rest of the methods see Ref. [57, 58]. 

There are two methods "under construction" which bear some promise to be suitable for direct 
QCD simulations at any \i ^ in the future. The first one is the method of stochastic quantization 
[59], where the evolution of the gauge fields is governed by complexified Langevin equations, 
which can be integrated even at nonzero \x. In practice the method is prone to instabilities and still 
at the stage where the conditions for the convergence to the correct ensemble are explored. This 
year results [^(J from testing the method in an SU (3) spin model, which approximates the strong- 
coupling large-quark-mass regime of QCD, show that the method gives the correct solution for this 
particular case. The application of stochastic quantization also was studied in the chiral random 
matrix theory at finite density [pT[|, where analytical and numerical solutions can be compared. A 
second method for direct simulation and nonzero density is the world-line approach [|2|], the basic 
idea of which is to integrate out explicitly the gauge fields instead of the fermion ones as is usually 
done. This way the sign problem could be alleviated []63|]. Unfortunately the gauge field integration 
cannot be performed explicitly due to the four field interactions, except in the strong coupling limit 
of QCD. In this limit after gauge integration the remaining degrees of freedom are mesons and 
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Figure 9: (Left) The CEP temperature estimates from the first max of c„ [|68[]. Solid line is the full model 
Te- (Right) The CEP location determined with different methods and compared to possible positions of the 
freezeout curve. Taylor expansion method (1): standard staggered Nf = 2, N t = 6, volumes up to V = 24 3 , 



m„/m p w 0.3 [^6|; (2): p4 Nf = 2 + 1, Nt = 4, volumes up to V = 24 \ m n /m p « 0.3 ||70|]. Reweighting 

(3) : staggered Nf = 2+ 1, Nt = 4, volumes up to V — 12 3 , m K /m p — physical [pip. Canonical approach 

(4) : /?4 iV/ = 2,N t = 4, volume V = 6 3 , m K /m p ps 0.7 [Q. Direct canonical simulation (5): clover A'/ = 3, 
= 4, V = 6 3 , w s = 700 - 800 MeV 



baryons i.e., color singlets. The "worm" algorithm Q64p is one best suited for this setup and this 
year a new "continuous time" version J65l ] of it was introduced which does not suffer from the sign 
problem and was applied to the study of the phase diagram in the strong coupling regime of QCD. 
Both of these methods (stochastic quantization and world-line approach) were reviewed previously 



in more detail [57] 



4.2 The Taylor expansion method 

The main idea of the Taylor expansion method is that the basic thermodynamic quantities are 
represented as series whose coefficients can be calculated on an ensemble with pL = where a 
traditional MC simulation is possible, which is the method's most important advantage. In other 
words, the expansion is around jU = and the expansion parameter is /J./T. For the pressure the 
expansion explicitly is: 



p(T,n) _ InZ _ f, m /M\" m 1 d"lnZ 

T 4 VT 3 L> C "V){ T ) ' C »^) n\VT^d{y:/TY 



(4.2) 

11=0 



where pi is the light quark chemical potential and we assume the heavy quark one is zero for 
simplicity. 

The expansion coefficients c n (T) are nonzero only if n is even due to the CP symmetry of the 
partition function. The expansion is roughly expected to converge when p is smaller than the first 
Matsubara frequency i.e., when p<nT. However, in the chiral limit the radius of convergence of 
the Taylor expansion will shrink to zero at the CEP. Even at nonzero quark masses and at small p, 
in the crossover region the statistical fluctuations in the Taylor coefficients grow large and make 
them more difficult to determine. The method has also other drawbacks which should be taken 
cautiously into account. For example, the number of terms in a given c n grows approximately as 
6" and there are large cancellations between them. The c„'s are estimated stochastically and the 
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number of noise vectors required to keep the noise under control grows exponentially with n and 



the volume Q57|]. As a whole the computational effort grows factorially with increasing order of 
the expansion, the highest reached for QCD calculation being n = 8 [|66|]. Higher orders are also 
sensitive to finite volume effects since essentially c n is an «-point correlation function. 

The position of the CEP in principle can be estimated using the Taylor expansion method. The 
quark number susceptibility, whose expansion is 

ft (r,M)_^(r )M )_f n{n _ l)Cn{T) (tLY-\ (4 . 3) 



t 2 dQi/T? ^ — jJl^Vr 



is expected to diverge at the CEP (unlike the pressure which stays continuous). The radius of 
convergence of the series for the pressure and the Xq(T,n) is the same and can be related to the 
coordinates of the CEP (jj, E , T E ) if the following condition is fulfilled. If in the complex plane the 
closest singularity of the expansion happens at real /J. then the radius of convergence r is related to 
the CEP coordinates as 

1/2 

(4.4) 



He/T e = lim 



n- 



C2n+2 



1/2 

= lim 

n— >oo 



c x 



Mathematically, for the singularity to lie on the real axis there should exist an order no such that 
for each n > no, all c n > 0. 

An estimate for Tg can be obtained independently of r from the temperature behavior of the 
coefficients c n . Close to the CEP a derivative of c n with respect to the temperature is equivalent to 



two derivatives with respect to \x, essentially meaning that dc n /dT ~ c n+ 2 Q68|]. This property can 
help in the determination of the temperature around which all coefficients with order higher than n 
are positive, by finding the temperature of the first maximum of c n . This gives consecutive estimates 
for Te as n is increased which are expected to converge from above to the CEP temperature. From 



estimates of r and Te the critical chemical potential jj, E can also be found through Eq. ( |4.4[ ). The 
most important question about the Taylor expansion method is: how many orders are needed for a 
reliable determination of the CEP coordinates (or other thermodynamic quantities)? The answer to 
this question through direct QCD simulations is difficult due to the large computational resources 
required to reach higher orders. Instead some insights into this question may come from model 
calculations for which these type of computations are much cheaper. 

The PQM model, whose degrees of freedom are Polyakov loops, quarks and mesons, has been 
shown to be a good approximation of QCD for 2+1 flavors and at /J. = when solved in a mean field 
approximation $7^. This year the same model has been explored (still in the mean field approxima- 
tion) in the case of ji. ^ [j58|]. The authors compare the mean-field full solution results with ones 
obtained from a Taylor expansion to 0(22). In Fig. || (right) shown is the result for the distance to 
the phase boundary from successive estimates of the convergence radius, which is determined both 
from the coefficients of the pressure expansion and from the quark number susceptibility one (see 
Eq. (fO])). Also shown are data points for the first pole in Pade approximations to the pressure and 
X q (T,^l), which is an alternative method for finding the distance to the phase boundary. From the 
figure the following conclusion can be made: the estimates r\ n _ 2 converge faster to the full solution 
than the r2„ ones, and the fastest convergence is achieved through Pade approximation of % q (T,ii). 
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Figure 10: (Left) The critical line as determined in Ref. [ [75| | and its position with respect to the freezeout 
curve. (Right) The isentropic energy density for different constant values of the ratio of the entropy S and 
baryon number Ng [fjjj]. 



It appears that at least O(10) in the Taylor expansion of % q (T,pL) is needed in order for the result 
to be within 10% of the full solution one. In Fig. |9] (left) shown are the data for the estimate for 
the temperature of the critical end point T CM , as a function of the order of the coefficient c n , from 
whose maximum it has been determined. From the figure it follows that the approach to the full 
solution Te is very slow with this method and even for O(20) the estimate is still about 10% off, 
which can lead to a significant systematic error in the determination of This study may point 
to a potential problem with slow convergence of the Taylor expansion for the true QCD case, to 
the extent the model is a relevant description for it. It is also possible that the slow convergence is 
a feature of the model rather than the real QCD case, meaning that a particular model's behavior 



may only be appropriate just as a guide. Still, another recently studied QCD model Q69Q (3d Ising 
model) estimates that at least 0(16) of the Taylor expansion is needed, in support of the claim that 
the inclusion of high orders may be necessary. 

4.3 The critical end point and the curvature of the critical line 

In this subsection I review some recent quantitative results for the features of the QCD phase 
diagram at nonzero chemical potential. For example, the location of the CEP has been the object 
of a number of studies using variety of methods [^6, 7(J 71, [72|, 73]. In Fig. ^ (right) a number of 



recent values for it are shown. There are two results obtained from the Taylor expansion method 
as described in the previous subsection, denoted by (1) and (2). One of the results (3) is from 
reweighting techniques and the rest are from canonical approaches (points 4 and 5). The results 
appear to cluster between n/T c = 0.5 — 1.0 and T /T c = 0.9 — 1.0, except one of them (4). Also 
shown in this figure are a host of possible freezeout curves which differ only on the assumed value 
of T c to which they are normalized. It appears that the cluster of CEP results lies close to the 
freezeout curves. (However, the results for the CEP location from the analytic continuation method 



[ 74 ] are in a qualitative disagreement with all of the above. They are interpreted to mean that the 



CEP may not even exist at physical quark masses.) 
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Another important feature of the diagram is the curvature of the critical line, and more inter- 
estingly, its comparison to the one of the experimentally determined freezeout curve. The critical 
line curvature K at jit = can be determined using the Taylor expansion method where to lowest 
order 



T c (ixl) = T c (l-K-^/T c 2 



dT c (}i 2 B ) 



with = 3jJ. being the baryon chemical potential. In a recent work 
through the following expression: 



d<$> 



50 
df 



T c ,fi B =0 



(4.5) 

H B =0 

1 K has been determined 



(4.6) 



where <p is an auxiliary observable, for example {yy) or % q , the derivatives of which are calculated 
at the critical temperature and on an ensemble at }i B = 0. This study with 2+1 flavors of stout 
fermions was done at N t = 6,8 and 10 with physical quark masses and the continuum limit of the 
data was taken. The results for the curvature are fc^/r 2 ) _ 0.0089(14) and = 0.0066(20) 

from data for the strange quark number susceptibility and the chiral condensate. In Fig. [l0| (left) 
the resulting critical lines (following from both observables) are shown at small jj, B assuming that 
the curvature does not depend on the temperature. The two values for fc are compatible within their 
errors and are also broadly compatible with results from other studies. For example, a p4 study 



with 2+1 flavors gave K = 0.0066(5) using 0(4) scaling analysis [76]. Results with two flavors also 
appear close to these values: K = 0.0059(1) from analytical continuation with staggered fermions 
and K 0.0078 from 0(4) scaling analysis with Wilson fermions [21]c. This closeness 



between the 2 and 2+1 flavor results is based on the expected similarity of the critical behavior in 
both cases. The curvature of the freezout curve at }Xb = is experimentally estimated to be k 0.02 
[HJ], which is a few times larger than the estimated values for the critical line one, meaning that, 
provided the curvature of the critical line does not change with T, the two curves will increasingly 
diverge as jj, B grows. 



4.4 The EOS at nonzero chemical potential 

The EOS at nonzero jU is important for HIC experiments when the center of mass energy is low, 
in which case the baryon density grows. More accurately, the relevant EOS is the isentropic one, 



where the ratio of the entropy and the baryon number (S/Nb) is constant. In Fig. [lOj (right) results 



for the energy density at different S/Nb from one of the most recent studies Q45Q are shown. The 
values for S/Ng are chosen such that they correspond to the ones achieved at different experiments: 
30, 45 and 300 are for AGS, SPS and RHIC, respectively. The calculation was done with the Taylor 
expansion method to 0(6), 2+1 flavors of asqtad fermions at = 4 and 6 and mi/m s = 0.1. As 
the temperature is changed the light and heavy quark chemical potentials were tuned such that the 
strange quark number density is zero (within statistical errors). The results show that for large 
S/Nb the EOS is difficult to distinguish numerically from the pLj s = (i.e., S/Nb = °°) case. At 
smaller S/Nb (such as those relevant for AGS and SPS) the energy density deviates from the zero 
chemical potential case at temperatures above the crossover one. 
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5. Other topics in nonzero temperature and density 

The list of subjects covered in this review is by no means exhaustive and the chosen topics are 
only some of the ones closely related to the HIC experiments or cosmological observations. Here I 
list several of the topics and recent results, which I did not cover in my talk, and their contributors. 

• QGP transport coefficients: J. Langelage 

• Chiral magnetic effect: A. Yamamoto, T. Ishikawa 

• Dirac operator spectrum: F. Pittler, T. Kovacs, H. Ohno, Z. Lin 

• t/(l) restoration: P. Hegde (plenary), G. Cossu 

• String tension: P. Bicudo 

• Large N c or Nf. M. Panero, K. Miura, M. Ogilvie, N. Yamamoto 

• SU (2) studies: P. Giudice 

• Strong-coupling limit of QCD: O. Philipsen, S. Lottini, K. Miura, A. Li, etc. 

6. Concluding remarks 

During the last year there has been solid progress in the study of QCD at nonzero temperature 
and density. One of the important advances is the resolution of the discrepancies in the values 
of T c determined with different versions of improved staggered fermions. Another is that more 
understanding was gained of the possible electromagnetic and finite volume effects on the T c (the 
latter were found in the SU (3) case only). It will be useful to extend such calculations to the fully 
dynamical case and also to determine these effects on other quantities such as the EOS. 

There are some developments in the calculations of the EOS itself. The first one is its deter- 
mination in the case of 2+1+1 dynamical flavors. The second one is that a new method for the 
calculation of the EOS was introduced and results for the entropy density at three values of the 
temperature in the pure gauge case were presented. However, in the 2+1 flavor case, the HotQCD 
and the WB results for the EOS at temperatures above the critical remain quite different and cur- 
rently the reason for this is unknown. The ongoing calculations with the HISQ action at N t = 12 
may be helpful in the resolution of this problem. 

The mapping of the QCD phase diagram at zero and nonzero density continued this year and 
new results for its properties were presented. 

There is some progress in the development of new methods for direct simulation of QCD at 
nonzero ;U, such as the stochastic quantization and the world-line approach; however, these methods 
remain yet to be applied in the true QCD case. 

Numerous studies of QCD-like theories strive to achieve at least a qualitative understanding of 
the behavior of the true QCD theory in certain limits (such as the strong-coupling one) and many 
new results became available this year. 

In conclusion, QCD at nonzero temperature and density is a very dynamic area of research 
which encompasses a large number of subjects related both to experimental and purely theoretical 
problems. This review cannot do justice to all of them for which I apologize. I want to thank 
everybody who sent me their results in advance and were available to answer my questions. I also 
want to thank the organizers of this conference for providing me with the exciting opportunity to 
present and discuss these topics. 
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